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Abstract 

We do a semiclassical analysis for two or three spins which are coupled antiferro- 
magnetically to each other. The semiclassical wave functions transform correctly under 
permutations of the spins if one takes into account the Wess-Zumino term present in 
the path integral for spins. The Wess-Zumino term here is a total derivative which 
has no effect on the energy spectrum. The semiclassical problem is related to that of 
anyons moving on a sphere with the statistics parameter 6 being 2irS for two spins and 
3itS for three spins. Finally, we present a novel way of deriving the semiclassical wave 
functions from the spin wave functions. 
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This Letter illustrates three different ideas using a simple quantum mechanical 
model. These ideas, which have been studied recently in other contexts, include a 
semiclassical treatment of quantum spins for large values of the spin S [ 1 - 5 ], a 
Wess-Zumino (WZ) term which is present in the path integral for spins [ 2, 6 ], and the 
possibility of fractional statistics on the sphere [ 7, 8 ] . We consider two or three spins 
coupled antiferromagnetically. Although a naive semiclassical analysis reproduces the 
low-lying energies and degeneracies correctly, there is a marked difference between the 
ways in which the naive semiclassical and spin wave functions transform under the 
exchange of any two spins. The complete symmetry group is the permutation group 
Sn of N spins in the quantum theory, and the braid group on the sphere in the 
naive semiclassical theory. We relate this difference in transformation properties to a 
WZ term which appears in the coherent state path integral Lagrangian for spins. This 
term treats each spin as a charged particle moving on a sphere which has a magnetic 
monopole of strength 4nS at its center. Once this term is taken into account, the 
modified semiclassical wave functions have the correct transformation properties under 
spin permutations. We can think about the semiclassical problem in terms of N anyons 
moving on a sphere with the statistics parameter 6 (which is defined modulo 2ir ) 
being 2ttS for N = 2 and 3irS for N = 3. Thus the naive semiclassical wave 
functions for three half-integer spins exhibit semionic statistics. At the end, we will 
verify our analysis by directly constructing the correct semiclassical wave functions from 
the spin wave functions . 

Consider first an antiferromagnetic Hamiltonian for two spins H = (S\ + 
S2) 2 = P ■ For any value of the spin S, the energies are 1(1 + 1) with a de- 
generacy 2/ + 1 where / = 0,1,2, ... , 2S . The three-j symbols have the symmetry 
[9] 
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where mi , m-2 and m denote eigenvalues of ( Si , ( S2 )3 and Is respectively. 
Hence, under the exchange P12 = Si <-> S2 , the wave functions transform by the 
phase ( P12 ) qu = ( —1 ) 2S + 1 where the subscript 'qu' denotes quantum . 

Now consider a semiclassical treatment of this problem. For S >> 1 , we introduce 

— * — * — * 

a vector <$> = ( Si — £2 ) / 25 . This satisfies 

T-$ = 
and ? = ! + --_ 

— * 

For low-lying excitations (i.e. / << S ), we see from (1) that is an unit vector. 
The naive semiclassical Lagrangian is 

1 - 2 

C = - <j> (2) 

with the constraint 4> 2 = 1 . Canonical quantization of (2) reproduces the above 
Hamiltonian except that — / 2 is now given by the Laplacian V 2 on a sphere. The 
semiclassical energies are therefore /(/ + 1) with degeneracy 21 + 1 , and the naive wave 
functions are the spherical harmonics Y^ m (a, (3) . (Here (a, (3) are the polar angles 

— * 

such that 4> = (sin a cos (3, sin a sin f), cos a ) . Namely, the direction of spin (or 
particle) 1 has the coordinates (a, (3) while particle 2 is at the point (tt — a, it + (3) 
where <a<n, 0<{3<2n). Semiclassically, / can be any non- negative integer. 
This spectrum agrees with the exact one for < / < 2S . It is interesting that the 
semiclassical energies and degeneracies are correct even if S is not much greater than 

— * — * 

one and / is not much less than S. Under the exchange of the two spins, — > — (j). 
The wave functions Y^ m then transform by the phase ( P12 ) n .sc = ( — 1 ) 1 where 
the subscript 'nsc' denotes naive semiclassical . 
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The difference of (—I) 25 between ( P12 ) qu and ( P12 )nsc can be explained 
as follows. The two-spin problem can be semiclassically thought of as two particles 
moving on a sphere with the Hamiltonian forcing them to lie at antipodal points for 
low energies. Under an exchange , the two particles will together trace out a closed 
curve which encloses a solid angle 2ir. Now, it is known that the path integral for spins 
contains a WZ term which makes each particle see a magnetic monopole of strength 
4irS at the center of the sphere [ 2, 6 ]. A particle which goes around a closed curve 
enclosing a solid angle Q picks up an Aharonov-Bohm phase exp (iQS) . Thus an 
exchange of the two particles produces a phase exp (i2nS) . 

To be explicit, the WZ term equals S $ (1 + cos a ) for particle 1 if we choose 
the Dirac string of the monopole to pass through the north pole a = . The WZ 
term for particle 2 is then S j3 (1 — cos a ) . The sum of the two is 

C wz = 2S/3 (3) 

which is a total derivative. The correct semiclassical wave functions are therefore 
exp (i2(3S) Yi^m (a, (3) . Note that since the phase factor 77 = exp (i2(3S) is single- 
valued (except at the points a = or tt when one of the particles lies on the 
Dirac string), the energy spectrum is unaffected by (3). But the new semiclassical wave 
functions do show the correct exchange phase (Pi 2 ) qu due to the factor 77 . 

To summarize the semiclassical picture, the particles behave like anyons on a 
sphere with the parameter 6 = 2nS . For anyons on a sphere, 9 is only allowed to 
have the 2{N — 1) values given by np/ (N — 1) , where p = 0, 1, 2, ... , 2A — 3 [ 8 ] . 
So we may hope to find a truly anyonic behavior (i.e. 6*^0 or tt ) if N = 3. We 
therefore turn to the more interesting problem of three spins. 

The Hamiltonian H = & + S2 + S3) 2 = P can be shown to have the low-lying 
spectrum /(/ + 1) with degeneracy (21 + l) 2 for < / < S . (For / > S , the 
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expression for degeneracy is different). Here / takes integer or half- integer values if S 
is an integer or half-integer respectively. For any /, let us consider the 21 + 1 states 
which have the eigenvalue h = I , and study their transformation properties under 
the permutation group S3 of the three spins. (Since the group operations of S3 
commute with the total spin operators / , the same transformation properties will 
hold for other values of I3 also). We have only examined small values of S and /, 
namely, < I < S < 2 . For integer S, we find that these 2/ + 1 states fall into / + 1 
irreducible representations (IR) of S3 consisting of / doublets and one singlet. Under 
the exchange of any two spins, the singlet picks up the phase ( P ) qu = ( — 1 ) s + 1 . 
(For example, the ground state (1 = 0) has a totally symmetric wave function if 
S is an even integer, and an antisymmetric wave function if S is an odd integer). For 
half- integer S, the 2/ + 1 states with I3 = I fall into / + 1/2 doublets. Under any 
exchange , the doublets always transform by a 2 x 2 matrix whose eigenvalues are 
± 1 . 

Some of the above statements can be understood using the two-spin results. To 

— * — * 

obtain a total spin / with three spins, the spin of two of them, say Si and S2 , must 
add up to values lying in the range S — /,..., S + I . (This explains the (2/ + l)-fold 
degeneracy for a given value of /3 ). Under the subgroup P\2 , the exchange phases 
are therefore ( — 1 ) 35 ~ 1 , ... , ( — 1 ) 3S + 1 . It then follows that for half-integer 
S, there are / + 1/2 states with P\ 2 = — 1 and / + 1/2 with Pi 2 = 1 , 
while for integer S, there are / states with P12 = ( — 1 ) s + 1 + 1 and / + 1 with 
P12 = ( — 1 ) s + 1 ■ Of course, it needs more work to derive the singlet and doublet 
structure under the full group S3 . 

Now we do a semiclassical analysis. For large S, we introduce two vectors [ 4 ] 
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— * — * 

^ Si — S2 

^ . (4) 

r Si + S2 — 2S3 
and 2 = - 

One can again derive identities similar to (1) which show that for S much greater 
than both / and 1 , (f>( = 2 = 1 and 0i • 02 = . We introduce a third 
vector 03 = </>i x 2 and define an SO (3) matrix R whose columns are given by 
the three vectors 0i , 2 and 03 respectively. The total spin / transforms R from 
the left. Thus [k , R] = — Ti R where i = 1,2,3 and the matrices Tj have the 
elements ( T t ) jk = - i e ijk ■ 

There is a two-to-one mapping from the manifold of SU(2), namely the hyper- 
sphere S 3 , to the matrices R . (This is made explicit in Eq. (13) below). We 
parametrize points on S 3 by an SU(2) matrix V_ = X4 + i a ■ x where a are the 
Pauli matrices and x\ + x\ + x\ + — 1 . For later use, we define 

z\ = x\ + i X2 = sin a e tb 

(5) 

and Z2 = X3 + i X4 = cos a e tc 

where < a < n/2 and < b, c < 2n . The symmetry group of the manifold 
S 3 is SO(4) = SU{2) x SU{2) / Z 2 . These two SU (2) 's are generated by matrices 
I and / ' which act on R (or V) from the left and right respectively [ 10 ]. (The 
two can be thought of as rotations about a set of axes which is either space-fixed or 
body-fixed). On R , a rotation acts as an orthogonal matrix O = exp (ien-T) while 
on V_ , the corresponding action is by an unitary matrix U_ = exp ( ie h-a /2 ) where 
h is an unit vector and e is the rotation angle. (Notice that we are using the same 
symbol I to denote both the quantum spin operators and the semiclassical angular 
momenta which act from the left. It will be clear from the context which one we mean). 
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The semiclassical Lagrangian can now be shown to be [ 4, 11 ] 

£ = I tr R T R = \ tr V* V (6) 
8 2 

where i? T = RT 1 and v} = V^ 1 . Canonical quantization of this yields 
H = — V 2 / 4 where V 2 is the Laplacian on S 3 . The normalization in Eq. 

— * 

(6) is fixed by considering small fluctuations near the identity. R = Is + i 2 e-T 
corresponds to V_ = I2 + i e- a which is near the north pole of S 3 , namely, X4 — 1 . 
(Here I n denotes the n x n identity matrix). Then (6) becomes C = e , the 
Hamiltonian is II e 2 /4 , and the Laplacian on S 3 is V e 2 = — II e 2 . 
On S 3 , the left operators / have the form [ 10 ] 

h = lua + - 491 - zsas) = + I) (7a) 

l + — Z2 d* — z\ 82 and /_ = z\ 82 — z\ d\ 
while the right operators / ' are 

/3' = i(, 1 a 1 -, 2 9 2 -, 1 *a 1 * + 495) = |(|-|) (n) 

1+' = Z2 d\ — z\ 82 and = z* — Z2 d\ 

where di denotes d / dzi . The Laplacian is V 2 = — 4 / 2 = — 4 / /2 . The harmonic 
functions on S 3 transform as the representation ( I, I ) under SU(2) x SU(2) . Here 
/ can be an integer or half-integer. The degeneracy is (21 + l) 2 . (For instance, the first 
five harmonics are 1 for / = , and z\ , Z2 , z\ and z\ for I = 1/2 ). The Laplacian 
acting on these harmonics gives — 41(1 + 1) . Under a 2n rotation of the spins about 
any axis, R — > R , V_ — > — V_ and the harmonics transform by ( — 1 ) 21 . So we 
have to choose / to be an integer (or half-integer) if S is an integer (or half-integer). 
Thus the semiclassical theory reproduces the correct spectrum for < / < S . 
The proper manifold for semiclassical quantization is therefore SO (3) for integer S and 
SU(2) for half-integer S as we might have expected for an odd number of spins. 



The 2/ + 1 semiclassical wave functions with / 3 = / are zf , z 2l ~ x z<i , ••• , z 21 . 
An exchange of any two spins is equivalent to a 180° rotation about the third spin 
which is a matrix acting on R or V_ from the right. (For V_ , one has to further 
specify whether the rotation is clockwise or anticlockwise as the corresponding matrices 
differ by a minus sign). Under any one of the three possible 180° rotations, we find 
that 



where P is an off-diagonal unitary matrix whose eigenvalues are ± i , not ±1 . For 



and V_ from the right by i a 2 or — i a 2 . So the matrix in (8) is P_ 12 = i cr 2 or 
- i a 2 ■ 

Hence, for half-integer /, the 21 + 1 functions given above fall into 1 + 1/2 doublets 
each transforming irreducibly under the braid group B3 , namely, 



Under any exchange , these doublets transform with eigenvalues ± i . For integer /, 
on the other hand, the 21 + 1 states fall into / doublets as in Eq. (9), but there is also a 
singlet given by z[ z\ . Under any exchange , the doublets transform with eigenvalues 
± 1 while the singlet picks up the phase ( P ) nsc = { —1) 1 . Thus in all cases, the 
eigenvalues of the (one- or two-dimensional) exchange matrices ( P_^ ) nsc differ from 
those of ( P^ ) qu by the phase exp (i3nS) . 

We can understand this by using the WZ term and a typical exchange shown in 
Fig. 1. The antiferromagnetic Hamiltonian forces the three particles to lie 120° apart 




(8) 



example, an exchange of Si and £2 corresponds to multiplying R from the right by 





(9) 
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on a sphere. Assume that S3 is fixed along the south pole (0, 0, — 1 ) and exchange Si 
and S2 by a 180° rotation. The resultant closed curve on the sphere encloses the solid 
angle Q. = 2n ( 1 + cos 60° ) = 3n . Hence the Aharonov-Bohm phase picked up is 
exp (i3nS) . Note that the sense of the exchange is important if S is an half-integer. In 
this Letter, however, we will not be precise about whether the argument of the phase 
is 3irS or - SttS . 

The symmetry group of three particles placed 120° apart on the sphere is the 
braid group B3 . A one- or two-dimensional IR of B3 can be obtained by taking an 
IR of the permutation group S3 and making it anyonic by multiplying the appropriate 
matrices by exp (id) for a clockwise exchange of any two particles. What we have 
shown above is that the anyon parameter 6 is equal to ± 3irS . 

We now study how these naive wave functions get modified if the WZ term is 
included in the Lagrangian . We first parametrize the three classical spins S n in terms 

— * 

of the orthonormal vectors (f> n given by 

4>3 = ( sin a cos (3, sin a sin (3, cos a ) 

4>i = ( — cos a cos (3, — cos a sin (3, sin a ) (10) 

and 4>2 = ( sin /?, — cos (3, ) 



Then 



Fi 1 r, ■ 1 27m . -* , 27m . -> ,„ „ . 

S n I S = sin ( 7 + — ) fa - cos ( 7 + — ) <f> 2 (11) 



for n = 1, 2, 3 . Now we find that the WZ term obtained by adding up the 
contributions from the three spins is again a total derivative. Indeed, 



£wz — S — 
at 



In ( ^ ~ 1 6XP ( ^ 7) - ) 1 (12) 

' V tan 3 (a/2) + % exp (- i3 7 ) / J v ' 

Next, we rewrite the matrix R in terms of the S s coordinates (5). This takes 
the form 

1 — 2 ( x\ + x\ ) 2 ( X1X2 + X3X4 ) 2 ( X3X1 — X2X4 ) 
R = I 2 ( 27x2 - X3X4 ) 1 - 2 ( x% + x\ ) 2 ( X2X3 + x\X4 ) I (13) 

2 ( X327 + X2X4 ) 2 ( X2X3 — X1X4 ) 1 — 2 ( x\ + x\ ) 
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A comparison then shows that a = 2 a, (3 = b + c and 7 = c — b . The 
WZ term is thus 

£wz = - iS±\n( i ~ H 3 ) (14) 
at \ z* + 1 z\ ' 

Hence the correct semiclassical wave functions are obtained by multiplying the naive 
wave functions by the phase factor 

v = ( i~ iz K ) s (15) 

(Our choice for the direction of the Dirac string and the expression for rj seem to break 
rotational invariance. However, one can define new angular momenta by the unitary 
transformation / — > rj I rj ~ ). 

We see that rj is single-valued on S s (unless zf = i z\ which corresponds to 
one of the particles being at the north pole of the sphere) . Hence the spectrum remains 
unchanged. But the phase factor restores the correct permutation properties to the 
wave functions . For instance, under the exchange Pi 2 , z\ — > and Z2 — > — z\ . 
Then r\ changes by exp ( — i3nS) as required. (See the discussion following Eq. (21) 
below. If Pi 2 is taken to transform z\ — > — Z2 and Z2 — >• z\ , then r\ changes 
by exp (i3irS) ). 

Finally, we turn to a derivation of the semiclassical wave functions from the spin 
wave functions for S = 1/2 . For the direction n = ( sin a cos /?, sin a sin /3, cos a ), 
the eigenvector of h ■ a /2 with eigenvalue 1/2 is given by 

, . _ / cos (a/2) exp{-i(3)\ , , 

1 ' P * ~ \ sin (a/2) J {W) 

The phase is chosen such that (16) is ill-defined only at the north pole. The ket \ ) 
for the two-spin problem is defined to be the tensor product 

I $ ) = I a, (3 ) \n - a, n + {3 ) (17) 
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where the first and second factors denote the vectors (16) for spins 1 and 2 respectively. 
The exact ground state with / = / 3 = has the wave function 

= Mi) •(!)-£(!)•© 

We then define the semiclassical wave function to be the amplitude ( 4> | -0o,o ) which 
turns out to be exp lo,o («, 0) / V% ■ Similarly, the three excited states with 
I = 1 have the wave functions ( <fi \ ipi,m ) = exp (i(3) Yi )Tn / VQ . (The Y^ m are 

— * 

normalised according to the measure d<p = sin a da d(3 / 4tt ) . 

For the three-spin problem, we define the ket | V_ ) as the tensor product 

| V ) = | a x , 0! ) <g> | a 2 , fo ) ® \ a 3 , fo ) (18) 

where the three kets | a n ,(3 n ) can be deduced from Eq. (11). To be explicit, the 
bra ( a n ,(3 n \ takes the form 



where u> = exp (i27r/3) . It is then clear that any 3-spin wave function | vp ) will have 
an amplitude (V_\ ip ) which is a polynomial in Zi and z* multiplied by the phase 
rj in Eq. (15). For instance, consider the two ground states with / = / 3 = 1/2 . 
Their wave functions 

- £ (o) ' (o) • (?) 

(20) 

transform under the exchange P\ 2 with phases — 1 and 1 respectively. One then 
finds that 

Vs 

(V\ipi) = —(z 2 +iz 1 )r) 

% <21) 

and ( V | ip2 ) = — ( zi + i 22 ) »7 
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If Pi 2 takes z\ — > ^2 and 22 — ► — zi , then 77 — > i 77 (that is, 
77 — > 77 exp (— i3irS) ). Hence ( F | Vi ) an d ( H I "02 ) transform correctly 
under P i2 . 

This construction can be generalized to any spin S . Let us work in a basis in 
which ( S n )s is diagonal. The key ingredient is the eigenvector | a, (3 ) (with 
eigenvalue S ) of the matrix h • S . The m th entry of the column | a, (3 ) is given 
by [9] 

i^^ = /(^(»l)' + "(-»l)'""-- w - w < 22 > 

where m takes the values S, S — 1, ... , — £ from top to bottom. This is well- 
defined everywhere except at the north pole. Then the two-spin ket in (17) produces 
wave functions of the form 

( $ I ) = N s ,i, m e t2S ? Y hm (23) 

where Ns,i, m is a normalization constant. To prove (23), we observe that in terms of 
the Euler angles (a, f3, 7) , the vectors exp (iS/3) | a, (3 ) and exp {iS(3) \ % — a, it + (3 ) 
are given by the first and last columns of the rotation matrix (a, (3, 7) respec- 

tively. These two columns are proportional to exp (—iSj) and exp (iSy) . The 
7-dependence cancels when we take a direct product of the two. Thus exp [i2S(3) \ (J) ) 
is given by a (2S + l) 2 -dimensional column of the direct product (a, (3, 0) <g) 

(a, f3, 0) . The transformation properties of this object under rotations and 
consequently the differential equation satisfied by it are then precisely the same as 
those of the spherical harmonics Y^ m (a, (3) . (If there had been a dependence on 7 , 
then the transformations and the differential equation would have been different since 
they would involve d/d^y [ 9 ] ). Thus, 

e" l2Sp < $ I l\ ^i, m > = I e- l2S(i < $ I ^i,m ) (24) 
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On the left hand side of (24), I = Si + S2 is a matrix which acts on | ip^ m ) while 
on the right hand side, / is a differential operator acting on exp (— i2S(3) { <fi \ ipi,m ) • 
(For example, /3 = — i d/d(3 ). Then Eq. (23) follows from the fact that | -0;,m ) 
has P = /(/ + 1) and Is = m . 

The constants Ns,i, m = Ns,i in (23) are independent of m due to rotational 
symmetry. We can choose N$,i to be real and positive. The completeness of the 
I tjj^m ) and the fact that ( <fi \ <fi ) = 1 means that 

2S 

(2/ + 1) N%,t = 1 

1=0 

Note that the (2S + l) 2 x (2S + l) 2 matrix 

M — J d$ \$) ($\ 

is not proportional to the identity. The orthonormality of the | V^m ) implies that the 
eigenvalues of M are given by the constants Ng t . An explicit formula for Ns,i can 
be derived by considering the state with Is = I . The corresponding three-j symbols 
are given by 

1/2 



& & I \ _ (_ l) S+l-m 

m I — m —I ) 



2l\ (2S-l)\ (S + m)l (S + l-m)\ 
I J (2S + 1 + 1)! (S -m)\ (S-1 + m)\ 



Using this one finds that 

e" t2S ? ( $ I ^1,1 ) - Y it i (a, (3) 

and the constant of proportionality determines 

(2S)l (2S)l 

s ' 1 (2S + 1 + 1)! {2S-l)\ { ' 

For the three-spin problem, the ket in Eq. (18) produces the wave functions 
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( Y_ | ?pl,m ) which are given by polynomials in Zi and z* times the phase in (15). 
From (16), (19) and (22), it is clear that the polynomial is even for integer S and odd 
for half-integer S . An analysis similar to the one for the two-spin case will show that 
the polynomials are now given by harmonic functions on S 3 . We will however omit 
the details here. 

It is worth noting that the semiclassical analysis in this Letter holds for values of 
spin as small as 1/2 . It would be interesting to extend these considerations to more 
than three spins [ 5 ] or even to a spin chain. We often know, by either numerical or 
exact methods, the way in which the quantum ground state and low-lying excitations of 
antiferromagnetic spin chains transform under discrete symmetries like parity (defined 
as the reflection of the chain about one site). One might ask how this is related to 
the symmetry properties of the semiclassical field theories which are typically some 
non-linear sigma models. 
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Figure Caption 



1. Three particles placed 120° apart on a sphere. Exchanging particles 1 and 2 
keeping particle 3 fixed traces out a closed curve which encloses the shaded area. 
The complement of this area covers the solid angle 3n . 
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